Universality and non-universality in behavior of self-repairing random networks 
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We numerically study one-parameter family of random single-cluster systems. A finite- 
concentration topological phase transition from the net-like to the tree-like phase (the latter is 
without a backbone) is present in all models of the class. Correlation radius index vb of the back- 
bone in the net-like phase; graph dimensions - d m i n of the tree- like phase, and D m i B of the backbone 
in the net-like phase appear to be universal within the accuracy of our calculations, while the 
backbone fractal dimension Db is not universal: it depends on the parameter of a model. 

PACS numbers: 



Properties of disordered networks, such as porous ma- 
terials, are of great interest from both practical and the- 
oretical point of view (see [lj]). The standard percola- 
tion models (see 0,0) can not adequately describe many 
of these systems, in particular the "single-cluster" ones, 
which are obligatory connected (e.g., because of mechan- 
ical stability requirements) and do not allow for detach- 
ment of finite clusters. In our earlier paper [3| we have 
introduced and studied a simple model of self-repairing 
bond-percolation (SRBP) in which all finite clusters aris- 
ing in the process of gradual destruction of randomly 
chosen bonds in the system are immediately "repaired" 
by means of regeneration of that very bond whose de- 
struction has lead to violation of connectivity. The sys- 
tem therefore remains connected and constitutes a single 
large cluster at all accessible concentrations. The SRBP- 
model turned out to be analytically tractable by means 
of a partial mapping onto the standard bond-percolation 
model; we were able to prove the existence of a finite- 
density topological phase transition, at which the back- 
bone of the system vanishes and the system occurs in a 
fragile tree-like state with zero mechanical rigidity and 
zero electric conductivity. The backbone density Pb and 
correlation radius £b slightly above the transition follow 
the scaling laws 

Pb(p) oc (p - Pc f* , Zb(p)<x(p-P c )- Vb . (1) 

The critical concentration of bonds p c and relevant criti- 
cal exponents could be related to those of standard bond 
percolation (though not being necessarily equal to them) . 
In particular, in two-dimensional lattices 

^(SRBP) _ ^(pcrc) = 4757(10); (2) 
^(SRBP) s ^ (pcrc ) _ 4 ^ ^ 

The numerical value for /3^ >C1C ' ) in 2D, given in©, is 
taken from Ml, while the exact value of 
was first obtained in Q. The fractal dimension of the 
backbone is generally related to (3b and vb 

Db = D-Pb/v b . (4) 



For the SRBP model it is also the same, as in percolation. 
In 2D: 

D (SRBP) _ ^(perc) = ± g^g) (g j 

From the topological point of view the difference be- 
tween the tree-like phase and the usual net-like one is the 
abundant presence of arbitrary large (up to the system 
size) cycles in the graph of the net-like system and ex- 
ponential decay of concentration of large cycles on the 
graph of the tree-like system. 

One can easily envisage many other models of single- 
cluster random networks. We were not able, however, to 
develop an analytical approach to any physically reason- 
able model except SRBP. The question about universal- 
ity of the discovered phase transition (and, moreover, the 
very existence of the transition in a particular model) is 
very important and nontrivial. In contrast to standard 
percolation models, the connectivity constraint (common 
for all single-cluster models) is essentially nonlocal. This 
increased complexity of the model makes the existence 
of relations of single-cluster models to local field theories 
(similar to the well-known relation between the standard 
percolation and the Potts-model, see || and Q) highly 
improbable; we cannot therefore proclaim the universal- 
ity basing on standard renormalization-group arguments. 
Note, that the invasion percolation (see ]1,[3|), which also 
involves a nonlocal constraint in its definition, does not 
belong to the standard percolation universality class and 
is characterized by its own critical indices. Thus any 
new single-cluster model introduces a new puzzle of its 
own, and it is very interesting and important to explore 
different patterns of their behavior and try to establish 
some order in the corresponding zoology. In this paper 
we start this program, and study the non-universality of 
critical exponents. 

Our first step was numerical exploration of the self- 
repairing site-percolation (SRSP) on a square lattice, 
where randomly chosen sites of the lattice (not bonds, as 
in SRBP) are gradually removed and regenerated each 
time when the removal causes disconnection of a finite 
cluster from the mainland. Numerical simulation reveals 
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a well-defined phase transition at the site-concentration 
x c = 0.608(1) with the critical index of the backbone 
density /3^ RSP ^ = 0.463(2), which is distinct from the 
result @ for the SRBP-model. The latter difference is 
a definite manifestation of the non-universality of the 
phase transition. On the other hand, our simulation 
did not show any reliable difference between the SRSP 



index of correlation radius: 

j/SRBP) _ 

we obtain 
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1.337(8) and 



4/3. For the fractal dimension, using 
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1.653(2), 



(6) 



also different from Dg RBP ^ . 

Another important fractal characteristic of the critical 
backbone on the spatial scale R <C C is the graph dimen- 
sion D m in, describing the dependence of average "chemi- 
cal distance" £(R) (that is, the length of the shortest path 
between two sites of a backbone, separated by euclidean 
distance R) 



i B (R)ozR I 



As we have shown in 



D (SRBP) _ D (perc) _ 
min min 



1.13(2), 



(7) 



(8) 



where the numerical value (for 2D) was taken from Q. 
Note that the graph dimension of the percolation back- 
bone coincides with that of the entire infinite cluster 
(see 0). Our simulations for the SRSP-modcl give 
-^min SP ' = 1-136(10), that, again, is not different from 
the SRBP value (|5J) within our accuracy. 

In contrast with the SRBP-model, where the tree-like 
phase can only exist in the rangeptrec < P < Pa bounded 
from below by ptrce > (see [I), in SRSP-model the 
quasi-tree cluster can have arbitrary low concentration 
(0 < x < x c ). For x — > the concentration of cy- 
cles vanish very rapidly, so that the patterns looks like a 
tree already slightly below x c . In FigfT] reasonably large 
randomly chosen fragments of a quasi-tree samples with 
x = 0.5 and x — 0.25 are shown. It demonstrates that 
the cycles die out very fast with lowering x: the sample 
with x = 0.25 is already practically a tree. 

The tree-like phase is fractal in the entire range x < x c , 
so that the chemical distance £t(R) in this phase obeys 



e T (R) oc R d - 



1 < d min < 2. 



(9) 



As was shown (both numerically and analytically) in [J], 
for SRBP-model d m i n does not depend on p throughout 
the tree-like phase and exactly coincides with the similar 
index rf| I ^ T '' of the Minimal Spanning Trees ensemble 
(MST) on a lattice: 



,(SRBP) 



,(MST) 





FIG. 1: SRSP-patterns in the tree-like phase: upper panel, 
x — 0.5, here few small cycles are still visible; lower panel, 
x — 0.25: concentration of cycles is already so low, that we 
do not see any of them in this fragment of a quasi-tree. 



The latter numerical value of d 



(MST) . 



in 2D was obtained 



1.22(1). 



(10) 



in [10|,[llJ (see, also [4|). 

To determine d m i n and D m in for SRSP model in this 
paper we simulated systems of size 1024 x 1024 with open 
boundary conditions. Statistics on average Euclidean dis- 
placement (R) as a function of chemical distance £ was 
gathered by performing breadth-first searches on points 
from the central part of the system. About 10 7 starting 
points were processed, and plain power fits (R) ~ £ 1 / dmin 
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FIG. 2: Log-log plots for root mean square euclidean displace- 
ment R vs. chemical distance I for the SRSP-model on square 
lattice with four different concentrations: x = 0.62 (above 
the threshold); x = 0.608 (at the threshold); x = 0.57 and 
x = 0.20 (both below the threshold). Within the accuracy 
of our calculations, the slope of the first curve corresponds 
to d m in = 1.0002 (practically, unity), the second curve - to 
D m ix, — 1.136(10). The third and the fourth curves both have 
the same slope, corresponding to d m ln = 1.226(12) (practi- 
cally, the same as d^^ T ^). 

were applied in the stable region 500 < I < 900 not 
significantly affected by the boundary. Our simulations 
give, again, x-independcnt 

d£f P) = 1.226(12), (11) 

throughout the tree-like range of x. This value, within 
our accuracy, does not differ from rf| 1 ^ T ' ) . To our knowl- 
edge, however, there is no obvious relation between qua- 
sitrees, appearing in the tree-like phase of the SRSP, and 
the MST ensemble. 

At very low concentration i«l the tree-like phase of 
the SR(S/B)P-modcls apparently acquire a new spatial 
scale £c(zO ~ x~ vc » 1, - the density correlation length, 
that diverges as x — > 0. On the scales 1 -C L <C £c 
the quasi-tree, besides being a "chemical fractal" (with 
nontrivial d m i„ > 1) becomes also a "density fractal" with 
nontrivial fractal dimension Dq (see FiglSJ) ■ Properties of 
this low-density phase and corresponding critical indices 
will be studied in a separate publication. 

For further study of the non-universality we have con- 
sidered a one-parameter family of hybridized SR(S/B)P- 
models, in which, at each step of the process, with prob- 
ability 1 — Q a randomly chosen bond is removed (and re- 
stored, if necessary) and, with probability Q, a randomly 
chosen site together with all adjacent bonds is removed 
(and also restored, if necessary). In Fig[4]the concentra- 
tion dependence of the backbone density Pb (p) is shown 
for five different values of the parameter Q. For three of 
them (Q = 0, Q = 0.5, and Q = 1) the critical exponents 
Pb-i v and u> were accurately extracted from the data. 

We have devised an efficient algorithm for simulations 




FIG. 3: SRSP-pattern at low density x = 0.0125: not only a 
tree, but apparently also a density fractal. 
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0.4757(10) 


1.6432(8) 


1.85(1) 


0.5 


0.420476(1) 


0.468(2) 


1.649(2) 


1.785(5) 
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0.373116(1) 


0.463(2) 


1.653(2) 


1.788(3) 



TABLE I: Nonuniversal characteristics of the family of 
SR(S/B)P models. 

of models with arbitrary Q based on dynamic maintain- 
ing of connectivity on the dual lattice that allowed site or 
bond deletion operations to be performed in 0(1) time. 
Lattices of size 128-256 turned out to be the most useful, 
and 10 6 -10 7 samples were simulated for each Q. This 
algorithm enabled us to accurately determine fine differ- 
ences in the critical exponents, appearing in the problem. 

The procedure of finding the indices was as follows (c.f. 
paper [12J, where a similar procedure was proposed for 
study of backbones for standard percolation): For a given 
L an ensemble of realizations on L x L square was gener- 
ated with an additional "boundary constraint" , requiring 
that the sites and bonds, belonging to two opposite sides 
of the square (say, upper and lower ones) could not be 
removed and constitute two "bars" which by convention 
belong the the backbone. For a given realization of the 
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TABLE II: Universal indices of the SR(S/B)P- family. 
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FIG. 4: The backbone density for SR(S/B)P hybrid models 
at different values of mixing parameter Q. Concentration p is 
defined as a number of bonds present in the system; both p 
and Pb are normalized by the total number of bonds in the 
full lattice. 



process the threshold p c was defined as a concentration, 
at which the infinite cluster splits into two parts - one, 
connected to the upper bar, the other - to the lower 
bar. A distribution function of p c over the ensemble was 
found. The average and the standard deviation of this 
distribution are (see, e.g., Q): 

%(L) =J c ^ Pc + CiL~", A Pc (L) « CsiT 1 /^, (12) 

and the average value of the backbone density at the 
transition is 



P B (L) = Pb(Pc) ~C 3 L 



D b -D 



(13) 



where C\, C2, C3 - critical amplitudes. The val- 
ues of Pc, vb, Db, Pb, and to, given in the tables, 
were extracted from relations (|12I13[) . The extrac- 
tion procedure is illustrated in FigO where /3 c s(L) = 
(v/\na) \n[P s{aL) / P b{L)\, and a > 1 is some rescaling 
parameter (the result is not sensitive to the choice of a). 
Note, that the exponent u) is not universal already within 
the class of standard percolation (namely, it depends on 
a particular lattice, as well as on the boundary condi- 
tions), so that the dependence of w on Q may only be 
viewed as an indirect evidence of non-universality of the 
SR(S/B)P class. 

In conclusion, we have studied a one parametric family 
of self-repairing site/bond percolation single-cluster mod- 
els. In all models of this family a topological phase tran- 
sition between net-like and tree-like phases was found, 
but the backbone fractal dimension Db{Q) turned out to 
be non-universal. It apparently is a smooth monotonous 
function of the site/bond mixing parameter Q, varying 
by ~ 0.01 in the range Q < Q < 1. Though the above 
variation is quite small, it is well outside the error bars 
±0.002 of our calculations. To reach this high accuracy 
special efforts were made. 



FIG. 5: Linear extrapolation of p e s(L) to 1/L 
rescaling parameter was chosen a = 3/2. 
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The non-universality of the critical indices was not un- 
expected, since there are no known reasons for universal- 
ity, similar to those, existing for the family of standard 
percolation problems. In this connotation, it was a great 
surprise that index vb of the backbone correlation ra- 
dius; the graph dimension D m i n of the backbone; and 
the graph dimension rf m j n of the tree-like phase seemed 
to be Q-independent within our (rather high) accuracy. 
Based on our numerical observations, we conjecture: 

-,(SR(S/B)P), 



(SR(S/B)P; 
J B 



(Q) 



(pore) 



D, 



,(SR(S/B)P) 



(Q,p) = D 
(Q,p) 



(pcrc) 
min ' 
,(MST) 



So far, we did not find any rational explanation for this 
intriguing phenomenon of " partial universality" . 
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